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The Pfaffian quantum Hall states, which can be viewed as involving pairing either of spin- 
polarized electrons or of composite fermions, are generalized by finding the exact ground states 
of certain Hamiltonians with k + 1-body interactions, for all integers k > 1. The remarkably 
simple wavefunctions of these states involve clusters of k particles, and are related to correlators of 
parafermion currents in two-dimensional conformal field theory. The k = 2 case is the Pfaffian. For 
k > 2, the quasiparticle excitations of these systems are expected to possess nonabelian statistics, 
like those of the Pfaffian. For k = 3, these ground states have large overlaps with the ground states 
of the (2-body) Coulomb-interaction Hamiltonian for electrons in the first excited Landau level at 
total filling factors v = 2 + 3/5, 2 + 2/5. 
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I. INTRODUCTION 

Trial wavefunctions for fluid states with all the parti- 
cles in a single Landau level (LL) have played a paradig- 
matic role in the fractional quantum Hall effect (FQHE) 
|jj since the work of Laughlin [Q. The Laughlin states 
were extended to other filling factors using the hierarchy 
approach which can also be used to generate trial 

wavefunctions, though these are not as elegant or unique 
as Laughlin's. A different motivation for extensions of the 
Laughlin wavefunctions to yield other filling factors is the 
composite fermion approach [pi, which generates wave- 
functions that are simple, but again not quite unique for 
each filling factor. In some cases, trial wavefunctions are 
also energy eigenstates of lowest energy for some Hamil- 
tonian. In the FQHE, these Hamiltonians have usually 
been found to annihilate the appropriate trial states, and 
whether this is so can be established without much diffi- 
culty Some early examples are in Refs. j3|||. 

However, these distinct trial wavefunctions at a given 
filling factor do not necessarily correspond to distinct 
phases of matter. Phases of matter should be character- 
ized by their ground state quantum numbers (including 
filling factor or Hall conductance) , types of long-range or- 
der, ground state degeneracy (if any), and properties of 
excitations such as whether or not there is an energy gap, 
and the quantum numbers of the elementary excitations 
that serve as building blocks for all others. In brief, these 
characteristics involve universal aspects of the asymptotic 
low-energy, long-distance physics, and not details of mi- 
croscopic wavefunctions, either trial ones or exact en- 
ergy eigenstates of a realistic Hamiltonian. One usually 
tries to argue that the model Hamiltonian for which the 



trial ground and some excited states are exact eigenstates 
yields properties that are characteristic of a phase. If it is 
true that the energy spectrum has a gap, and hence that 
the fluids in question are incompressible, then the uni- 
versal properties should be insensitive to perturbations 
of the Hamiltonian, and the properties found from the 
trial states do typify a phase. However, the existence of 
a gap is hard to establish for the model Hamiltonians, 
since it involves eigenstates of nonzero energy, which are 
not known exactly. For the Laughlin state and the pseu- 
dopotential Hamiltonian of Haldane ||, the arguments 
are fairly convincing. (We note that there are also com- 
pressible liquids |?],|| , which we believe also to be distinct 
phases, but for which there are so far no Hamiltonians for 
which exact eigenstates are known. The same comment 
applies to crystalline phases.) 

In the lowest (or M = 0) Landau level (LLL) in 
single-layer systems, the picture that has emerged is that 
the incompressible fluids are in the phases typified by 
the Laughlin and hierarchy states, or by the compos- 
ite fermion approach, these two approaches yielding the 
same phase for each filling factor with an odd denom- 
inator ^,0. (There are complications involving spin 
pl| , which we will not discuss here.) In higher LLs, the 
pscudopotentials j^], that characterize the interactions 
within the highest partially-filled LL, arc different, and 
the ground states must be reconsidered (the relative im- 
portance of the Zeeman energy may also change, but will 
not be considered in this paper). For example, the quan- 
tized Hall plateau (i.e. incompressible fluid) observed at 
i/ = 5/2 = 2|l/2@, where the LLL is filled with elec- 
trons of both spin components, and the 1/2 is the filling 
of the first excited (or N = 1) LL, does not correspond 
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to a hierarchy state; no corresponding plateau is seen at 
v = 1/2 in the LLL. A trial state in which the spins of 
the electrons in the J\f = 1 LL were unpolarized was pro- 
posed 1^, however, recent work suggests that the true 
ground state may be polarized ]l^ , p^ |. Moreover, even 
for simple fractions such as 2 + 1/3, the ground state for 
the Coulomb-interaction values of the pscudopotcntials 
in the W = 1 LL is on the borderline between compress- 
ible and incompressible states, and the Laughlin state 
does not have such a good overlap with the exact ground 
state, even when the V\ pseudopotential has been in- 
creased so as to enter the incompressible region ||. 

It is clear that the study of other trial wavefunctions, 
that are hopefully representative of distinct phases from 
the hierarchy, could still be of further use in under- 
standing incompressible fluids in higher LLs. One class 
of these is the paired states, in which one attempts to 
form Laughlin states of pairs of electrons, in some sense 
^5|,|6| Jl^Jl7[| . In particular, the so-called Pfaffian state 
was introduced and related to correlation functions in 
two-dimensional conformal field theory (CFT) by Moore 
and Read |ll| . This state occurs for filling of the topmost 
LL of v — l/g, for q even for fermions, the physical case. 
Using this correspondence, they argued that this incom- 
pressible fluid phase has fractionally-charged excitations 
of charge l/2q in electron units, and that these possess 
nonabelian statistics, a generalization of fractional statis- 
tics that we will discuss later. Also, there are neutral 
fcrmion excitations; the fluid ground state can be viewed 
as a Bardeen-Cooper-Schrieffer paired state of these neu- 
tral or composite fermions, and the fermion excitations 
are created by breaking pairs. The Pfaffian ground state 
is an exact, zero-energy eigenstate of a certain Hamil- 
tonian containing a three-body interaction p7| , p8| . The 
claims about the statistics were reinforced by later work 
which exhibited the gapless fermion excitations at the 
edge [ ^9|j20| l , and the degeneracy of the quasihole states, 
for the three-body Hamiltonian |H],[l8|. In the recent 
work |l3|]l4f |, the true ground state for the Coulomb in- 
teraction in the M = 1 LL was found to have a sizable 
overlap with the Pfaffian state which can be increased to 
large values (97 %) as V\ or V3 are varied about their 
Coulombic values. (Morf jl3| also found a large overlap 
with another paired state.) Therefore it is possible that 
the state observed at v = 5/2 is in the phase described 
by the Pfaffian, as first suggested in Ref. jl7j . 

In this paper we obtain a class of trial wavefunctions 
by a direct generalization of ideas that are valid for the 
Pfaffian state, namely the zero-energy eigenstates of a 
k + 1-body 6- function interaction, where k is an integer. 
The existence of the states is established by arguments 
using operator product expansions (ope's) in CFT p2[ , 
and the explicit form of the wavefunctions is obtained and 
proved to give the ground states. The states involve clus- 
ters of k particles, generalizing the pairs in the Pfaffian 
state. Further analysis gives predictions about the quasi- 



hole states and nonabelian statistics, which we partially 
confirm by solving the special Hamiltonians numerically. 
Then, in Section III we analyse the ground state of the 
Coulomb interaction in the Af = 1 LL for v = 2+3/5, and 
compare it with our state at the same filling factor; ex- 
cellent agreement is found, indicating that our functions 
are serious contenders to describe the phases in higher 
LLs, at least in some cases. Further discussion is given 
in the Conclusion. Some of our results appeared in an 
earlier unpublished work [ p3[ . 

II. SOLUTION OF SPECIAL K + 1-BODY 
INTERACTION HAMILTONIANS 

In this Section, we show how to solve certain interac- 
tion Hamiltonians, in the sense of finding the zero-energy 
eigenstates, with the help of operator product expansions 
(ope's) in conformal field theory (CFT) @. We first 
reinterpret the Pfaffian state in this light, then discuss 
our Hamiltonians. Then we analyse the quasihole states, 
and discuss nonabelian statistics and the number of sec- 
tors of edge states. 

A. Notation and Hamiltonians 

We will first define some notation for a system of par- 
ticles in the lowest Landau level (LLL) on the sphere jl| 
(it was used previously in Ref. Q). The magnetic field 
is radial and uniform with a total of flux through 
the surface, and in the lowest Landau level (LLL) each 
particle has orbital angular momentum N$/2. The LLL 
wave functions on a sphere are usually written (in a cer- 
tain gauge j|) in terms of "spinor" (or "homogeneous") 
coordinates Ui and Vi for each particle i = 1, . . . , N, 
with Ul = e 1 ^/ 2 cosOi/2, Vi = e^/ 2 sin 0j/2, where 
9i, are the spherical polar coordinates on the sphere. 
Since these imply that it;, Vi are not independent com- 
plex numbers, it is often more convenient, and will sim- 
plify the writing, to use a nonredundant parametriza- 
tion of the sphere by a single complex variable. This is 
done by stereographic projection, which gives the defini- 
tion Zi — 2Rvi/ui, where R is the radius of the sphere. 
Single-particle basis states in the LLL then take the form 
zf I '(1 + N 2 / '4i? 2 ) 1+Ar * /2 > wh ere the L z angular momen- 
tum quantum number is L z = N^/2 — m, so to < N^. 
Many-particle states can thus be written as 

* = *[](l + |z J | 2 /4i? 2 )-( 1 +^/ 2 ) (2.1) 

i 

and W must be a polynomial of degree no higher than 
N$ in each z^. Therefore, in the following we need spec- 
ify only 5» in order to describe a state. The function 
^ for a ground state on the sphere can also be used 
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to construct a wavefunction suitable for a corresponding 
disk-shaped ground state on the plane, by multiplying by 
exp(-i£,N 2 ). 

The Pfaffian state is defined as |l6| 

* Pf (zi, ...,z N )=Pi ( —^—) TT(«i - zj) 9 . (2.2) 
The Pfaffian 

Pf M v - = (2 jy / 2 (iV/2)!)- 1 £ sgnP [J M P(2r _ 1)P(2r) 

(2.3) 



i<3 



N/2 



P r=l 



of an antisymmetric N x N matrix M (N even) is the 
antisymmetrized sum over all pairings (zi — Zj)~ x (the 
analogous pairing in the spin-singlet case appears in the 
Haldane-Rezayi state ||). The filling factor of the state 
isv = N/N^ -> 1/q as iV -> oo (since A^ = q(N-l)-l), 
and g is odd for a boson state and even for fermions. 

The Pfaffian state for q = 1, q = 2 is the ground state 
of a three-body Hamiltonian • For these cases, the 

Hamiltonian penalises the closest approach of any three 
particles allowed by the statistics. Thus for q = 1, where 
the particles are bosons (note that we reserve the term 
"particles" for the underlying particles, which are either 
charged bosons or charged fermions [electrons], and not 
for composite particles), the Hamiltonian can be taken 
to be flfl 



H = V 



E 

i<j<k 



S 2 (zi - Zj)6 2 (zi - zk), 



(2.4) 



where the sum is over distinct triples of particles [Q . For 
numerical purposes on the sphere, it is more convenient 
to work with a projection operator form of the three- 
body Hamiltonian, instead of the 5-functions in (2.4). 
The closest approach of three particles on the sphere cor- 
responds to the state of maximum possible total angular 
momentum for the three. If the particles are bosons, the 
largest possible total angular momentum is 3A^/2 (recall 
that each particle has angular momentum N^/2). Then, 
for the q — I case, the Hamiltonian may be taken as pro- 
portional to the projection operator onto the (unique) 
multiplet of maximum angular momentum for each triple 
of bosons: 



H= Y, VP ijk (3Nt/2). 



(2.5) 



i<j<k 



The same trick works for the three-body interaction of 
fermions giving the q = 2 case; in this case, the maximum 
total angular momentum of three particles is 3N ( p/2 — 3. 
On the plane, the latter Hamiltonian corresponds to 
derivatives of delta functions. For these two cases, the 
Pfaffian state on the sphere is the unique energy eigen- 
state of zero energy at the stated N<j, value. (We will 



refer to such states as zero-energy states hereafter.) For 
larger q, these Hamiltonians can be generalized, in such a 
way that the zero-energy states are obtained from those 
for q = 1 by multiplying by Y[( z i ~ z j) q ~ 1 (it is assumed 
that for q odd, we are discussing bosons, and for q even, 
fermions). The presence of the latter factor implies that 
they are all zero-energy eigenstates of the projection op- 
erators for any two particles onto relative angular mo- 
mentum M = 0, 2, . . . , q — 3 (q odd), or M = 1, 3, 
. . . , q — 3 (q even) [or the corresponding total angular 
momenta N^, — 2, . . . , — q + 3, (q odd), etc.]. 
The space of states annihilated by such projections is in 
one-to-one correspondence with the full space of states 
of the q = 1 case, and the desired three-body projection 
operator [onto angular momentum — 3(q — I)] is 

the unique one that corresponds under this mapping to 
that already mentioned for q = 1. For each q, the Hamil- 
tonian can then be taken to be the sum of the three-body 
and all of these two-body projection operators. A very 
similar approach works for the other Hamiltonians stud- 
ied in this paper, so that results for higher v~ x can be 
deduced easily from those for the minimal v^ 1 for each 
type of state. These Hamiltonians can also be written in 
terms of 5-functions and their derivatives, so as to arrive 
at a form suitable for use in geometries other than the 
sphere. 

The goal of this paper can now be stated: we wish 
to generalise the solution for the zero-energy states of 
the three-body Hamiltonian (2.4) or (2.5) to k + 1-body 



Hamiltonians of the same closest approach form, for ev- 
ery k (note that for k — 1, we can consider the solutions 
to be the Laughlin ground and quasihole states ||). We 
will demonstrate the existence of an elegant solution to 
this problem by using conformal field theory, construct 
the ground state wavefunctions explicitly, and then com- 
pare the states with solutions for the Coulomb interaction 
in the second LL. 



B. Ope's and the Pfaffian state 

First we will reconsider the Pfaffian state and show 
how its property of being the unique zero-energy state 
for a certain Hamiltonian, that is of vanishing as three 
particles come to the same point, is related to operator 
product expansions in the corresponding CFT. We then 
generalize this to find the solutions to the simplest exam- 
ples of k + I-body Hamiltonians, for each k. The use of 
CFT leads to an existence proof for these wavefunctions, 
and in principle determines the wavefunctions uniquely. 

The basic idea of Ref. ]l6| is that the wavefunction of 
our 2 space-, 1 time-dimensional system can be related 
to a correlator of a certain conformal field theory. Here 
we are not interested in the simple Laughlin- Jastrow fac- 
tors, but in the other parts that may produce nonabelian 
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statistics. For the Pfaffian case, this part is the theory of 
free massless Majorana Fermi fields in 2 Euclidean space- 
time dimensions lira: 



* Pf = (^( Zl )...^))J|( 2 



Zj)\ 



(2.6) 



in which the correlator can be evaluated using Wick's 
theorem and (ip(z)ip(w)) = — (ip(w)ip(z)} = {z — w) ■ 
The Fermi fields can also be characterized by the operator 
product expansions (ope's) which hold inside correlators 



ip{z)I{w) 



(z - w)- l [I + 2(z - w) 2 T{w) + 
ip(w) + [z — w)dip(w) + . . . , 



(2.7) 
(2.8) 



as z — > w, where / is the identity (which is a constant), 
and T(z) = —^ipdip is the stress tensor. Throughout, 
d(f>(z) will mean dcf>(z)/dz and dots . . . will denote higher 
order, less-singular terms. While the correlator of Majo- 
rana fermions can be found by Wick's theorem, the use of 
ope's is more general and can be applied in a vast number 
of situations where there is no simple Wick's theorem. 

We now show how the ope for the fermions guarantees 
that the state is a zero-energy eigenstate for the three- 
body interaction of Greiter, Wen and Wilczek |L7| . This 



can be seen directly 17 from the explicit wavefunction 
( |2.2| ), however, the ope's provide a general argument that 
can be used even when the wavefunctions are unknown, 
if the function is a correlator of fields with known ope's. 
Suppose that in the correlator (without loss of gener- 
ality) first zi — > zi, then 23 — * z\, and take the most 
singular term of each product. The first limit (using Eq. 
fl2.7p) gives {z\ — z 2 )~ 1 I 1 and the second (using Eq. (fES)) 
then gives (z% — z\)~ 1 ip(zi). Multiplying by the Laughlin- 
Jastrow factor 'J'lj = Ili^ ^i — z j) we find that for 
q = 1 vanishes as (Z3 — Zi)(z3 — z 2 ), i.e. quadratically. 
Hence for the short-range interaction among bosons in 
Eq. (|]J) we obtain zero energy (ItJ . Furthermore it is 
the densest such state [ ^7| , ^o|Jl^1 . The same argument 
applies for q > 1 by construction. 



C. Parafermion states 

We next generalize the ideas to solve Hamiltonians 
with k + 1-body interactions. The simplest such Hamil- 
tonian, which will give the highest density such state for 
each k, which is always a state for bosons, is a (S-function 
interaction between the k + 1 particles: 



H = V 



i\<i2<—,ik+i 



S 2 (z n - Zi 2 )8 2 {zi. 



(2.9) 



As we will see, the wavefunctions involve dividing the particles into clusters of k each. To generalize the approach just 
described for the Pfaffian, we will consider a natural generalization of the ope algebra of Majorana fermions, which is 
the algebra of parafermions p5[. This consists of a set of fields rpt(z), £ = 1, . . . , k— 1 with ope's (ip e (z) = ipk-i(z)) 



M*)M*') - c ^'( z ~ z')- {Ae+Ae '- Ae+e,) l^+fAz') + ■■■} {£ + £' <k), 

M^W) - ct, k -t>(z - z')-^+^'- A '-^[^. t {z') + ...} (£'< £), 



2A f 



~( z ~ zT 2Ae [I +—(z- z'fT{z') + ...], 



T(z)Mz') 



At 



(z ~ z>) 



j^M*') + — f tol>t{z') + 



(2.10) 
(2.11) 

(2.12) 
(2.13) 



The consistency of the algebra puts conditions on the conformal weights Ag of the ipg's. The simplest solution to the 
conditions is Ag = £(k — £)/k, and with this choice the values of the central charge c = 2(k — l)/(fc + 2) and the 
numerical coefficients C£^< are then uniquely determined by consistency. With these choices, we obtain the algebra 
usually known as parafermions. (Other solutions for the A| are listed in Ref. p5[ , and some are studied further. 
In these cases the other parameters are not uniquely determined. These other algebras will not be needed here, but 
may be relevant in other FQHE problems. The procedure below will generate wavefunctions in all cases, but these 
may not always be solutions to simple Hamiltonians.) For k = 2, tj)\ — ip is a Majorana fermion, and the ope's are 
the same as Eqs. (2.7,2.8) above. 



Consider the function 



para \ 



■ ,z N ) = (ipi(zi) ■ 



I ( \\,T, M + 2 / fe 

"0lOiv))* L j 



(2.14) 



where N is divisible by k and M > is an integer. The ^^j +2 ^ k factor renders (2.14) nonsingular 

any Zi — ► Zj, 



(Zi 



according to Eq. (2.10) with A e + A& - A e+ £> = 2££'/k, which gives 2/k for £ = £' = 1. ^Ppara is totally 



symmetric for M even, antisymmetric for M odd, so describes bosons or fermions, respectively. As any Zj 



the correlator 



-2A, 



00, 



because the separations among the other N — 1 Vi's are then relatively small, and the ope's 
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imply that they "fuse" to form a single ip\. 
v = k/(Mk + 2) is the filling factor. 

Taking now the M = case, we can show that the function vanishes quadratically as any k + 1 particles come 
to the same point, say Z2, ■ ■ • , Zk+i approach z\ one by one, by a similar argument based on the ope's as for the 
Pfaffian. As the first k — 1 particles approach z\ , the most singular (or most slowly vanishing) terms must give the 
same result as the N = k case of "fpara, which is a constant (2f|; one additional particle gives ■>pi(z k+ i)I ~ ipi(zi), 
and the LJ factors give {zk+x — z\) 2 . Hence like the Pfaffian state with q = M + 1 = 1, this state is a zero-energy 
state for the k + 1-body <5-function interaction of bosons, and this fact again extends to M > 0. This shows that the 
desired wavefunctions exist, provided the CFT correlators do; but the parafermion theories are well-studied and can 
be constructed from other well-understood theories |25|]2"(| , so correlators with the properties specified by the ope's 
do exist. Notice that, had we chosen to use a different parafermion operator ipg (£ < k/2) in place of all the ipi's, we 
would have obtained a polynomial that vanishes more rapidly when k + 1 Zi's coincide, and a lower filling factor. 

Explicit functions can be obtained using the idea of grouping particles into clusters of k Jll| . The following procedure 
gives a symmetric *Ppara of the correct degree which is the Laughlin state for k = 1, and the Pfaffian state for k = 2: 
Divide the particles into clusters of k. For each pair of distinct groups, say z%, . . . , Zk and Zk+i, • ■ ■ , Z2k, we define 
factors x by (for k > 2) 



Hence ^p ar a is a polynomial of degree = v 1 A r — (M + 2), where 



Zfc+l, . . . , Z 2k ) = {Zl - Z k+1 ){z 1 - Z k+2 )(z 2 ~ Z k+2 ){z 2 ~ Z k+3 ) ■ ■ ■ (Z k - Z 2k )(z k - Zfe+l), 



(2.15) 



that is, each member of a cluster is connected with only two members of the other cluster, by a factor z% 



For 



k = 1, we would have \ — ( z i 
whole expression, to obtain 



=2 



Now we multiply these for all distinct pairs of clusters, and symmetrize the 



x para / _j 



]J X(zp(kr+1 
P€S N 0<r<s<N/k 



))•••, Zp(k(r+1)): Zp( ks+ 1), . . . , Zp(k(s+l))) 



(2.16) 



where Sjv is the permutation group on N objects, and 
the prime on the sum denotes that the summation can 
be restricted to permutations obeying P(l) < P(k) < 
■ ■ ■ < P(N — k + 1), which eliminates redundant permu- 
tations of the clusters; otherwise each term would appear 
(N/k)l times (in the expression for the Pfaffian above, 
we instead summed over all permutations but divided 
by the number of times each distinct product appeared) . 
Hence this function is totally symmetric, even with this 
restriction. To obtain larger M > 0, we must multiply 
by ^lj- The degree of the function in each Zi is then 
iV = M (N - 1) + 2(N/k - 1) = v- x N - (M + 2), as for 
the parafermion correlator times LJ factor above. 

Next we show that our function (for M — 0) vanishes 
quadratically whenever any k + 1 particles come to the 
same point. It is sufficient to consider each term in the 
sum over permutations separately. Clearly, we must then 
consider many possibilities, in which the k + 1 particles 
are distributed among different clusters, ranging from all 
but one in the same cluster, to each in distinct clusters. 
We will organise our proof in the following way: first we 
show that each term vanishes if k particles come to the 
same point, except in the case of all k in the same cluster. 
Then we use the evident fact that for k particles in one 
cluster, one in another, the term vanishes (this follows 
directly from the definition of \ itself). 

We must consider k distinct Zi's, chosen in any way 
from up to k distinct clusters. We may label these by the 
corresponding position in a cluster (from 1 to k), and the 



label (like r, s) for the cluster, that appear inside the per- 
mutation operator P. These determine when the product 
of x's vanishes. Put another way, it is sufficient to con- 
sider the term where P is the identity. Now consider the 
values for the position within a cluster as arranged on a 
circle (or clock face) , with the numbers increasing clock- 
wise until k is reached, which is followed by 1, so that 1 
is adjacent to 2 and k, 2 is adjacent to 1 and 3, and so 
on. The x's containing any two selected Zj's vanish only 
in two cases: (i) if members of distinct clusters occur at 
the same place on the clock; (ii) if members of distinct 
clusters occur at adjacent positions, but only if the posi- 
tion of the later cluster is arrived at by moving clockwise 
from the earlier by one step (see the definition of \) ■ We 
can then describe the structure of the term containing 
any set of k z^s by labelling positions on the clock face 
with the numbers of the clusters in which the Zj's ap- 
pear. Members of the same cluster must be at different 
clock positions because they are distinct ZiS. Clearly, if a 
term is to have any possibility of not vanishing, we must 
choose all the clock positions to be distinct. We then 
have a single number of a cluster assigned to each clock 
position. We will now show that it is not possible to ar- 
range these in such a way that the term does not vanish, 
except by choosing all the Zi's from a single cluster. We 
can view the clock as divided into regions (possibly con- 
sisting of positions that are not all adjacent) that have 
been labelled (or "colored") with a single cluster number. 
To avoid a vanishing factor, the numbers of the clusters 
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must decrease as one moves clockwise from one region to 
another. But the clock is not simply connected, and so 
this cannot be done all the way around the clock, unless 
there is only a single region (or cluster) involved. In the 
latter case, consideration of any one other Zi shows that 
the term vanishes, as mentioned before. Therefore, we 
have shown that the function vanishes at least linearly, 
but since it is totally symmetric after summing over per- 
mutations, it will actually vanish quadratically as any 
k + 1 particles come to the same point. 

We have shown that this vanishing property follows 
from the parafermion ope's, which should determine the 
function, however the vanishing property alone may be 
only a necessary and not a sufficient condition to de- 
termine that the function is the parafermion correlator 
times ^lj- Nonetheless, since we have demonstrated it 
for an explicit function, and since this function seems 
to be unique, on the basis at least of numerical verifica- 
tion for k — 3, 4 (discussed further below), we conjecture 
that the parafermion correlators for all k are given by this 
construction (by dividing lp|a by ^^J k )- It should be 
possible to prove or disprove this statement by studying 
the properties of our function as various combinations of 
Zj's coincide, and comparing these with the ope's. This 
was done for various paired states in Ref. [27| ], but we 
will not consider it further. 



D. Quasihole states and nonabelian statistics 

In this Subsection, we argue that states containing 
quasiholes (carrying charge l/(Mk + 2) each) for the 
parafermion states with k > 2 can be constructed in anal- 
ogy with the PfafHan case, and that nonabelian statistics 
are expected in all cases. The arguments go as follows. 
We will again begin with the Pfaffian state. 

The two-quasihole wavefunction proposed in Ref. fltj ] 
was 

*Pf+2qh(^l, ■ ■ ■]W 1 ,W 2 ) = Pf " J W£j 

(2.17) 

where f(z i ,z j ;w 1 ,w 2 ) = {z i -w 1 )(z j -w 2 ) + (z i -w 2 )(z :j ~ 
Wi). It can be interpreted as the insertion of two spin 
fields <t(w): 

*Pf+2qh oc (tp(zi) ■ ■ ■ ip(z N )cr(wi)a(w 2 )} 

x *£j - ^i) 1/2 (^ ~ w 2 ) 1/2 . (2.18) 

i 

The spin fields induce square-root branch singularities in 
the fermi fields, described by the ope's J2^] 

ip(z)a(w) ~ (z - w)- 1/2 a(w) + .... (2.19) 



The branch singularities are cancelled by the explicit 
square roots, to ensure that the wavefunctions are single- 
valued. This fixes the charge of each quasihole to be l/2q. 
We note jl6jthat quasiholes can be created only in pairs. 

In Ref. ]16[ it was proposed to extend this by defin- 
ing wavefunctions for 2n quasiholes by inserting 2n spin 
fields. Since the ope's, by definition, describe short dis- 
tance properties of correlators independently of what 
other fields are present, our argument above implies that 
the quasihole wavefunctions so obtained will all be zero- 
energy eigenstates of the appropriate H% (this can again 
be seen explicitly in the n — 1 case above) . This is anal- 
ogous to the Laughlin states, where for a suitable pseu- 
dopotential Hamiltonian the ground state and states with 
quasiholes added are zero-energy states j|. However, for 
the present case, there is not just a unique state for each 
set of quasihole positions. Instead there are 2" _1 linearly- 
independent functions (conformal blocks) of the z^s for 
In spin fields. This number follows from the ope of the 
spin fields §|: 

a(z)a(w) ~( z - w)~ 1/s I + const(z - w) 3/s i)(w) + .... 

(2.20) 

Here the higher-order terms fall into two families since 
the behaviour of each term as z — > w differs from that 
of one of the first two terms by an integer power of 
{z — w). Then in a correlator containing a{w\)a{w 2 ) • • •, 
as w\ —> w 2 , u>3 —f W4,, etc we can choose a member of 
either set of terms for each pair of cr's, which would give 
2™ terms, except that the total number of r/>'s in the cor- 
relator must be even, so we get 2™ _1 blocks, for ./V (the 
number of ip insertions) odd or even. In principle, the 
ope's can be used to fix the actual form of the correla- 
tors for 2n > 4. In practice, it was more convenient to 
find all the zero-energy states by explicit construction, 
and it was shown that for fixed positions of the quasi- 
holes, the number of states is 2"" 1 Thus the 
results for the Pfaffian are in complete agreement with 
these predictions based on CFT ope's. 

These degenerate spaces of quasihole states are the ba- 
sis for the nonabelian statistics properties of the quasi- 
holes (and similar properties are expected for quasielec- 
trons, or for combinations of quasiholes and quasielec- 
trons). As the locations of the spin fields are ex- 
changed by analytic continuation, the conformal blocks 
are mapped to linear combinations of each other (mon- 
odromy) p2fl . The conjecture of Ref. |l6| was that when 
the quasiholes are exchanged adiabatically as in Ref. p8| , 
these functions exhibit nonabelian statistics, meaning 
that the effect on members of the space is a linear trans- 
formation, which can be described by a 2™~ 1 x 2 n_1 
matrix by choosing a basis, and that this is the same 
as the matrix obtained by analytic continuation. This 
replaces the usual Berry phase that describes ordinary 
(abelian) fractional statistics; matrices representing dis- 



G 



tinct exchanges do not usually commute, hence the term 
"nonabelian" . The existence of many blocks is thus a nec- 
essary but not sufficient condition for nonabelian statis- 
tics (the adiabatic exchange has not yet been calculated 
explicitly) @. 

We now discuss the extension of these results to the 
parafermion states for k > 2. In place of the spin fields 
for the Majorana fermion, the parafermion system has 
"fields" (chiral vertex operators) & m p5| , p6| , where 1 = 0, 
1, . . . , k, while m is a periodic variable with period 2k, 
so m = 0, 1, . . . , 2k — 1, and further / — m = (mod 2). 
In the ope's of & m , the m's add mod 2k, however, there 



k-l 
k+vn 



= $L are 



are also identifications which imply that $ 
the same operator, so the values of m can be restricted 
so as to get each of the k(k + l)/2 distinct fields once; a 
convenient way to do this is by restricting —l<m<l. 
In this notation, <&[] = = I and ^e-k = ®2t = ^l- 
The special cases 07 = <&\, I = 1, . . . , k, (o^ = I), are 
called primary fields for the parafermion algebra. a\ in a 
sense generates the whole set by repeated operator prod- 
ucts, and is a natural analogue for a, to which it reduces 
for k — 2. We propose to insert a number nk of ax's into 
the correlator of parafermion currents to obtain the basic 
quasihole states; these are zero-energy states for our spe- 
cial Hamiltonians by the preceding argument. We use o\ 
because this leads to the minimal charge for the quasi- 
holes. Thus our proposal for quasihole states is to use 



(M) , , 
para+qh( Z l>---; W l:---> W «fc) = 

(ipl(zi) ■ ■ -^l(2jv)ffl(wi) ' ■■<Tl('W n k)) 

N nk 

M + 2/ kTTTT( „._ Wp) X/k 



(2.21) 



i=i p=i 



In these states, N, 



4 = u- 1 N - (M + 2)+n. Once again, 
this will in general define a whole set of states (confor- 
mal blocks), not a unique state. In the last factor, the 
exponent 1/fc is chosen to cancel the branch singularity 
as any ipi approaches any o\ , which is determined by the 
ope, 



^(z)*K0) ~ *- 1/fc 3j(0) 



(2.22) 



(see Refs. |2^,g6|]). This is the weakest singularity for any 
of the & m , so gives the smallest charge for a quasihole. 
The charge on the quasiholes is determined entirely by 
the Laughlin-like part of the wavefunction, not the con- 
formal block of the parafermions, so the usual plasma ar- 
gument shows that the charge (number of particles miss- 
ing from the vicinity of w p ) is v/k = l/(Mk + 2) (use of 
the other spin fields ai in general give charges vljk, and 
by including additional Laughlin quasihole factors we can 
add integer multiples of v to these charges). Mk+2 is the 
denominator of v — k/(Mk + 2), however k and Mk + 2 
have a common factor if and only if k is even, and this 
factor is just 2 (e.g. for the Pfaffian, k = 2). Therefore, 



when k is even, the charge is fractionalized compared 
with the usual value in a Laughlin or hierarchy state for 
a spin-polarized single component system at the same 
filling factor, which is always 1/q for filling factor p/q, 
where p, q have no common factors. The next nontriv- 
ial example is k = 4, where with M = 1 for fcrmions 
we obtain a v = 2/3 state, or v = 1/3 by particle- hole 
inversion, with charge ±1/6 excitations. 

The explicit function for n = 1, that is k quasiholes, is 
the same as ^para, except that a factor of the form 

T(zi, ...,z k ;wi,.. .,w k ) = _Q(^i - Wi) (2.23) 

i 

for each cluster of k particl es is inserted inside the sum 
on permutations P in Eq. ( 2.16 ); the permutations act 
only on the z,-'s, not the w p . This generalizes the func- 
tion / above for the Pfaffian, except that once again the 
symmetrization is now done all at once by the sum over 
P, which can still be restricted as before. The resulting 
function is also symmetric in the w p 's. 

For more quasiholes, we do not have the explicit func- 
tions in general, but we can count the number of zero- 
energy states for fixed positions of the quasiholes in the 
above construction, using the CFT, in a similar way as 
for the Pfaffian state above; we will consider k = 3 ex- 
plicitly. We require the ope's of the fields & m , which are 
given in Ref. |2l| for all k, and will not be written down 
explicitly here. Using the ope's one finds by repeatedly 
taking ope's with a±, that the number of conformal blocks 
in the parafermion theory for k = 3 for 3n quasiholes is 
a Fibonacci number, p3 n _2 ; where we define -Fi = 1, 
F 2 = 2, F 3 = 3, Fi = 5, and F m = F m ^ + F m _ 2 in 
general p0| . Thus, the number of linearly- independent 
zero-energy states for our Hamiltonian is also (at least) 
F$n-2 f° r fixed positions of 3n quasiholes, in the k = 3 
case, provided N is sufficiently large. For large n, this 
number approaches ~ (2 + \/5) n . As k increases, the ex- 
pressions for the n dependence for the parafermion cases 
will become progressively more complex than that for the 
Pfaffian (k = 2) case, which was 2 n ~ 1 . 

The results for the quasiholes can be compared with 
numerical results for the four-body Hamiltonian. The re- 
sults for the total number of zero-energy states (compare 
Ref. pq] for results on the Pfaffian state), were calculated 
for N — 6, 9, and 12 fermions at v = 3/5, but should be 
independent of M. First, the ground state at the stated 
is unique, and for small sizes was verified to be given 
by the explicit polynomial above. For one flux added 
(n = 1), that is 3 quasiholes, the number of states was 

N/i + 3 

as one would expect for 3 bosons in N /3 + 1 orbitals, or 
from the explicit wavefunctions above. For n = 2 flux 
added, that is 6 quasiholes, the numbers of states were 
in agreement with the formula 
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^N/3 + 6 



N/3 + 5 
6 



(N/3 + 4 
V 6 



(2.24) 



which is similar in form to those found for the Pfaffiaii 
state [ JL8| . Assuming this works for all larger N, and 
dividing by the first term which is the value of the po- 
sitional degeneracy that would be expected if the quasi- 
holes obeyed abelian statistics jl^] , we obtain the ratio 5 
in the thermodynamic limit, which agrees with the num- 
ber expected from conformal field theory. Similarly, for 
9 quasiholes (n = 3), the results agree with the formula 



N/3 + 9 
9 
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N/3 + 8 
9 



10 



N/3 + 7 
9 



(2.25) 



and so we expect the ratio to be 21 in the thermodynamic 
limit, as expected for fixed positions of the quasiholes. 
For more quasiholes, because of size limitations we have 
not been able to obtain any such formulas. Thus we find 
a satisfying agreement with our prediction, which tends 
to confirm that all the zero-energy quasihole states are 
obtained by inserting cti's in the parafermion correlator. 

We have not obtained results for the 4-body Hamil- 
tonian ground states on the torus, or the edge excita- 
tions (compare pb"||I§| ]), either analytically or numeri- 
cally. However, we expect that these calculations would 
lead in general to the conclusion that the number of sec- 
tors of edge states, or ground states on the torus, is 
(k + l)(Mk + 2)/2 (notice that this integer is divisible 
by the denominator of the filling factor in all cases, as 
required by Ref. pi]]). This result is based on a natural 
structure for the CFT, including the U(l) charge sector 
as well as the parafermion sector. Our analysis of this 
theory, into which we will not go in detail here, also in- 
dicates that while for k = 2, there are neutral fermion 
excitations, both at non-zero energy in the bulk and as 
gapless excitations at the edge |TqJl7|^o| , for k > 2 the 
analogous parafermion excitations carry charge 2v/k 
plus multiples of v. The excitation containing a tpi and 
charge 1 is identified with the physical hole, as usual, and 
in the CFT interpretation that applies (for example) to 
the CFT of edge excitations, it generates the chiral alge- 
bra, as in the Pfaffian case [ [ll3|j20j| . There are, however, 
neutral parafermion excitations for k even, that contain 
ipk/2- These can be viewed as being made from k/2 par- 
ticles and (Mk + 2) /2 flux. These are simply fermions 
for k/2 odd, bosons for k/2 even, and are like the usual 
composite particles. Also, for all k, there are nontriv- 
ial neutral excitations originating from the spin fields, 
= ®k m t ne earlier notation. For k — 3, the lat- 
ter, <&g, is the only field other than the parafermions ipe 
and the spin fields o - ;. These have no analogues either 
in the hierarchy or the Pfaffian (k = 2) states, though 
they do in the generalized hierarchy states, such as the 
331 state |2fl]. We further note some isomorphisms of the 
algebraic structures. For M = 0, the full chiral algebra, 
including the U(l) charge sector, is the level k SU(2) 



Kac-Moody current algebra, with representations (sec- 
tors) labelled by "spin" j = 1/2 = 0, 1/2, . . . , k/2 (this 
spin is of course not the physical spin, which is always 
polarized) . For M = 1, the chiral algebra is the so-called 
N = 2 superconformal algebra in the antiperiodic sector, 
and the representations for each k make up the known 
discrete series k = 1,2, . . . , for this algebra p^ |. In these 
cases, we have simply recovered known constructions of 
these algebras and representations from the parafermions 
p5| , p3[ . The special cases k = 1 (the Laughlin state) and 
k = 2 (the Pfaffian state) of these were mentioned earlier 
and the SU(2) k = 2, M = case was used in a 
recent paper [34]. 



III. COULOMB INTERACTION IN THE FIRST 
EXCITED LANDAU LEVEL 

We next turn to finite-size calculations. We have nu- 
merically constructed the wavefunctions of parafermion 
states for v = 3/5 (k = 3, M = 1), and v = 2/3 (k = 4, 
M = 1) for small sizes on the sphere and have confirmed 
that the k + 1-body Hamiltonians possess unique zero- 
energy ground states at the given flux N$. For k = 3, we 
have also obtained the excitation spectrum, both for the 
model k + 1-body interaction and the Coulomb potential 
in the J\f = 1 Landau level, and have studied the over- 
laps of the Coulomb ground state with our state as well 
as with the usual hierarchy states as the pseudopoten- 
tial V\ is varied about the its (AT = 1) Coulombic value. 
Below we describe these results. 
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FIG. 1. The low- lying spectrum for the four-body Hamil- 
tonian for TV = 15 electrons. The ground state is at L = 0, 
AE = 0. 

Fig. [l] shows the low-lying spectrum for N — 15 and 
v = 3/5 for the four-body Hamiltonian. The ground state 
is at L — 0, zero energy. The low- lying spectrum bears 
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some resemblance to the "hanging chain" shape seen in 
paired systems. 




FIG. 2. Pair correlation function of the v — 3/5 
parafermion state for TV = 18 electrons on the sphere, r is 
the great-circle distance. 

In Fig. |^ we show the pair correlation function for this 
same state but for the N — 18 size system, plotted as 
a LLL wavefunction, which describes the correlations of 
the guiding-center coordinates of the particles. Again, 
the "shoulder" at small r appears to be a characteris- 
tic of states in which the particles form clusters, such as 
pairs, and is also present for the Pfafnan state |^8| , p5[ . 
In addition, the large-distance oscillations are strongly 
damped. This feature indicates that the state is incom- 
pressible. 

We have also compared these states with the ground 
state obtained by diagonalizing the Coulomb potential 
for the TV = 1 Landau level. In the LLL, the hierar- 
chy states are most stable and these new ones will not 
be competitive. The first excited Landau level may be 
quite a different matter since V\ is reduced compared to 
V 3 (see e.g. Ref. M). Some typical values for Vi, V3, and 
V 5 at 7V = 22 are Af = 0: 0.4681, 0.2998, 0.2422, and 
J\f = 1: 0.4716, 0.3711, 0.2800. For the N = 1 Coulomb 
interaction ground state, the hierarchy state has small 
overlap-squared, whereas parafermionic states have very 
large ones. (The hierarchy state was obtained as the 
ground state of a model pseudopotential consisting only 
of non-zero V\ . We caution that these two states occur 
at different values of N$ because of the finite shifts on 
the sphere [JV^ = 5N/3 + 1 for the hierarchy].) For the 
k = 3 parafermion state we find 97% for N = 15 (where 
there are 36 states in the L = Hilbert space) and 88% 
for N = 18 (319 L = states), compared to at best 1 or 
2 % for the hierarchy 3/5 state. This is noteworthy since 
for these sizes we are very close to a Laughlin 2/3 state 
(i.e. the Laughlin 1/3 state of holes), because of the finite 



shifts in the N-N<p relations on the sphere. For N = 15 
the N<p for our state coincides with that of a single quasi- 
particle excitation of the v — 2/3 fluid while for N = 18 
it is at the same flux as the 2/3 condensate itself. This 
is a clear disadvantage for the parafermion states partic- 
ularly if we try to vary the short-range component of the 
Haldane pseudopotential V\ . Not surprisingly, increasing 
V\ by a few percent seems to favor the Laughlin state. 
However, one would expect the hierarchy v = 3/5 state, 
against which our state will be ultimately competing for 
large sizes, to show a slower rate of stabilization upon 
increasing V± away from its second Landau level value. 

To study this issue we compare the overlap-squared of 
the ground state of H = i?Coul,AT=i + SVi with both our 
state and the hierarchy 3/5 state. Fig. || shows these 
overlaps as a function of SVi for the N — 12 (52 L = 
states) hierarchy and N = 15 parafermion state. (The 
sharp drop in the parafermion curve at 6V1 about 0.03 is 
due to a level crossing: for larger 8V\, the ground state 
has L = 2.) It appears that our state remains stable for 
increases of 8V\ of up to 7-8% from the Coulomb value. 
Note however that the hierarchy state is not fully stabi- 
lized until V\ is increased by 20% of its Coulomb value. 
For large sizes there will be less interference from v = 2/3 
and the stability domain for our state may extend well 
beyond 7-8%. 
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FIG. 3. Overlap-squared of the two reference states, the 
v = 3/5 parafermion state (N = 15) and the v = 3/5 hierar- 
chy state (N = 12), with the state obtained by diagonalizing 
the M = 1 Coulomb potential with an added 8V1 component. 

Finally, in Fig. ^ we show the low- lying spectrum for 
the pure Coulomb case for N = 15. Again one finds 
some similarity with that of the four-body Hamiltonian 
(Fig. [I]), although in neither case is there a clear gap to 
a continuum of excited states. We defer precise gap es- 
timations incorporating finite layer thickness and other 
effects that modify the gap values, as well as detailed 
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studies of the quasiparticles, to future work. 
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tals and the usual hierarchy states (if we consider vary- 
ing the short-range part of the interaction away from its 
physical value, at fixed v). As the LL index Af increases, 
larger clusters with k ~ AT are expected to be favored 
[37|. It is very interesting that recent experimental work 
[38 1 has observed highly-resistive, highly-anisotropic be- 
havior, with nonlinear current- voltage characteristics, at 
total filling factors greater than 4, that is in the range 
Af = 2 to about 6. This behavior, seen around the center 
of each LL for each spin component, appears consistent in 
principle with the notion of a uniaxial or "striped" crys- 
tal phase, which was also predicted in Ref. |37| in this 
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region (the triaxial crystals of clusters were predicted for 
v away from 1/2), and thus may support these physical 
pictures. There is clearly much still to be done to under- 
stand the physics in this regime, in which we hope that 
the parafcrmion liquid states will play a role. 



FIG. 4. Same as Fig. [I] but for the Af = 1 Coulomb po- 
tential. The ground-state energy has been subtracted in the 
spectrum. 



IV. CONCLUSION 

To conclude, we have obtained the ground-state wave- 
functions for the k + 1-body Hamiltonians, and for the 
case k = 3, we have counted quasihole states, and (for 
M = 1) found excellent overlaps with the ground state 
of the Coulomb interaction for spin-polarized electrons 
in the first excited (Af = 1) LL, at total filling factor 
v = 2 + 3/5; by a particle- hole transformation, this also 
applies at v = 2 + 2/5. 

Our parafermion states contain clusters of k particles. 
For the Pfaffian state (k = 2) at v = 1/2, it has been sug- 
gested HQ that it may be favored in the Af = 1 LL, be- 
cause of this feature. That is, the correlation hole around 
each particle that is obtained by the Laughlin-style corre- 
lations, as a result of attaching q vortices to each particle 
(for filling factor of the topmost LL equal to 1/q) |36| , 
may not be sufficient to obtain the lowest energy in higher 
LLs, because of the "form factor" associated with these 
LLs — the real space wavefunction does not vanish when 
particles coincide G3j. However, by forming a cluster of 
two or more particles, the energy gained from the larger 
correlation hole surrounding the cluster may outweigh 
that lost within the cluster itself. Very similar physics 
is suggested by recent work motivated by the form of 
the pseudopotential interaction in higher LLs; the latter 
are also a consequence of the form factors. These works 
consider the formation of crystalline phases in which the 
particles are clustered, so there is more than one per 
unit cell p7| . If this idea is correct, then liquids contain- 
ing clusters might also be expected to occur in higher 
LLs, possibly as intermediate phases between the crys- 
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